AKARI NEAR-INFRARED SPECTROSCOPY OF 
SDSS-SELECTED BLUE EARLY-TYPE GALAXIES 
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ABSTRACT 

A near-infrared (NIR; 2.5 — 4.5/«m) spectroscopic survey of SDSS(Sloan Digital Sky Survey)- 
selected blue early-type galaxies (BEGs) has been conducted using the AKARI. The NIR spectra 
of 36 BEGs are secured, which are well balanced in their star-formation(SF)/Seyfert /LINER 
type composition. For high signal-to-noise ratio, we stack the BEG spectra all and in bins of 
several properties: color, specific star formation rate and optically-determined spectral type. 
We estimate the NIR continuum slope and the equivalent width of 3.29 /xm PAH emission. In 
the comparison between the estimated NIR spectral features of the BEGs and those of model 
galaxies, the BEGs seem to be old-SSP (Simple Stellar Population)-dominated metal-rich galaxies 
with moderate dust attenuation. The dust attenuation in the BEGs may originate from recent 
star formation or AGN activity and the BEGs have a clear feature of PAH emission, the evidence 
of current SF. BEGs show NIR features different from those of ULIRGs, from which we do not 
find any clear relationship between BEGs and ULIRGs. We find that Seyfert BEGs have more 
active SF than LINER BEGs, in spite of the fact that Seyferts show stronger AGN activity than 
LINERs. One possible scenario satisfying both our results and the AGN feedback is that SF, 
Seyfert and LINER BEGs form an evolutionary sequence: SF — > Seyfert — > LINER. 



Subject headings: galaxies: elliptical and lenticular, cD 
galaxies: active 

1. INTRODUCTION 

Early-type galaxies are often regarded as the 
objects at the final stage of galaxy evolution. 
Throughout many observations, it has been re- 
vealed that most nearby early-type galaxies have 
red colors and very poor contents of cold gas 
'e.gj" 

m- 



— galaxies: evolution — galaxies: formation — 

ditional star formation in the future. That is, most 
early-type galaxies seem to evolve passively, unlike 
late-type galaxies with active star formation. 

(|l999h 



However, since I Abraham et al 



and 
Menanteau et al.l (Il999t ) found a considerable frac- 



. iFaber fc Gallagherl [l976l iTrager et al.ll2000t 
Treu et alj|2006t . Such observational evidence ir 
dicates that early-type galaxies are composed of 
old stars and that they probably will not have ad- 



tion of early-type galaxies with blue colors and 
evidence of current star-formation in the Hubble 
Deep Field ( Williams et al.lll996 ). the existence 
of blue early-type galaxies (hereafter, BEGs) has 
become a new important factor to be considered 
in the formation scenario of early-type galaxies. 



BEGs are known to have positive color gradi- 
ents (i.e. center bluer than outskirt) on average 
and consist of both old stars and young stars 
( Menanteau et al. 2001at Elmegreen et al. 2005; 



Ferreras et alj|2005t iLee et al.ll2006U2008ll2010al) 
It has been also revealed that the internal color 
distributions o f BEGs are neither homologous 



nor s ymmetric ([Menanteau et al.l l2004t ILee et al. 



2006(), and that some BEGs are suspected to host 



activ e galactic nuclei (AGNs) ([Menanteau et al 



l2005t ILee et al.l 1200(1 l2008h . Unlike typical red 
early-type galaxies (REGs) preferring high density 
environments, BEGs tend to reside in intermediate 



or low density environ ments (jKannappan et al 



20091: iLeeet al.l IJOlObh . The origin and future 
of BEGs are still controversial. Observational 
evidence indicates that some BEGs may origi- 
nate from recen t galaxy interaction or mergers 
([Lee et al.1120061 ). but the gas inf all is also a possi- 



ble mechanism to build BEGs (Menant eau et al 



2001bl ) . After some time passes, BEG s ma 



evolv e into typical REGs (|Lee et al. 2006, 200 



ay 
)7, 



2008), but it is also possible that they evolve into 



bulges of late-type gala xies ([Hammer et al.ll2005r 
Kannappan et al.l 120091 ) . In any case, BEGs are 



probably in the forming/growing phase of early- 
type galaxies or bulges. 

To understand better the formation and evolu- 
tion of BEGs, we need to inspect what happens 
in the inside of BEGs. Since more than a half of 
BEGs have evidence o f their merger/interaction 
origin (JLee et al.l 120061) , the interaction- induced 
star formation (SF) is expected to make BEGs 
blue. On the other hand, galaxy mer ger /interaction 
some times also causes AG N activity (|Kewley fc Dopita 
120031; ISanchez et al.ll2005t IComerford et al.ll2009h . 
and the power-law co ntinua of AGNs m ay affect 
the color of BEGs. ILee et al.l (|2008l) reported 



that non-passi ve BEGs in the S loan Digital Sky 
Survey (SDSS; lYork et alJlJOOoh include both SF 
galaxies and AGN host galaxies, the number ra- 
tio of which is about 3:5 for -22.9 < Q1 M r 
111 < —20.7. However, optical spectra are eas- 
ily obscured by dust, causing the missing of a 
significant fractio n of AGNs in optical surveys 
([Best et al.ll2005l ). Moreover, ground-based spec- 
troscopy (e.g. the SDSS spectroscopy) of bright 



1 Absolute Petrosian magnitude in the r band with K - 
correction as if the object were at z = 0.1 iLee et al.l l 2008). 
All magnitudes in this paper are in the AB system. 



galaxies in the optical band often suffers the dilu- 
tion of low-co ntrast emission lines due to strong 
emission lines ( Ho et al.lll997l ). which also causes 
inaccuracy in estimating AGN activity. 

One of the most powerful tools to probe the 
SF/AGN activity is the spectroscopy in the in- 
frared (IR) b and. Polycyclic aromatic hydrocar - 
bons (PAHs; Is^iigren1ll984l: ILeger fe Pugetlll984l ) 
have been detected in various objects, such as 
Hn region, proto stellar clouds, planetary nebulae , 
and SF galaxies (|Peeters et alj|2002l; lHony||2002l) . 
It is known that UV photons excite the vibra- 
tional and stretching modes of PAH molecules, 
from the relaxation of which broad IR emission 
lines (such as 3.3, 6.2, 7.7, 8.6 and 11 . 2 jum 
lines) are produ ced (jCherchneff fc Barkerl Il989t 
Shan et al.l Il99ll ) . Such PAH emission features 



are go od indicato rs of SF activity. On the other 
hand, IVoiti (|l992l ) argued that PAHs in the cen- 
tral region of an AGN can be destructed by the 
strong X-ray emission from the AGN. Further- 
more, AGN activ ity is known to suppress SF itself 
(AGN feedback: lAntonuccio-Delogu fc Silkl 12008 ; 
Raffertv et al.l 120081 ). Thus, those properties of 



PAH emission are believed to provide a reason- 
ably clean diagnostic between SF and AGN. In 
addition, the shape of the IR continuum is an- 
other indicator of SF/AGN, because SF galaxies 
and AG Ns have different spec t ral energy distribu- 
tion (e.g. Ilmanishi et al.ll20 08; Risa liti et al.l l2006. 
2010l ). However, no systematic IR spectroscopic 



surveys of BEGs are seen in the literature. 

In this paper, we report the result of the first 
BEG spectroscopic survey in the near-IR (NIR) 
band, using the AKAR I (iMurakami et ail 120071 ) 
InfraRed Camera (IRC; lOnaka et al.l l2007l). The 
outline of this paper is as follows. Section 2 de- 
scribes the target selection and the AKARI spec- 
troscopic observation. Section 3 explains the re- 
duction and analysis of the AKARI data. The ex- 
tracted NIR spectra of BEGs and their features 
are shown in §4. The results are discussed in 
§5, and the main results and their implication 
are summarized in §6. Throughout this paper, 
we adopt the cosmological parameters: h = 0.7, 
0\ = 0.7, and fl M = 0.3. 
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Fig. 1. — Basic optical-NIR properties of the ob- 
served 36 blue early- type galaxies (BEGs). (a) 
Absolute magnitude in the r band as a function of 
redshift; (b) Apparent magnitude in the K s band 
as a function of redshift; (c) the color - magni- 
tude relation; and (d) the color - color relation. 
The small dots are SDSS early-type galaxies and 
the large symbols are the BEGs: stars for star- 
forming (SF) BEGs, triangles for Seyfert BEGs 
and circles for LINER BEGs (classified in Fig.©. 



Fig. 2.— The distribution of SDSS early-type 
galaxies (small dots) and the BEGs (filled circle) 
in the line flux ratio diagram, (a) The line is 
the boundary dividing SF galaxi es (lower-left) and 
AGNs (upper-right), given by iKauffmann et al. 
( 20031 ). (b) The line is the boundary dividing 
Seyferts (upper-left) an d LINERs (lower-right), 
given by iKewlev et al.l (|200 6). The number of 
filled circles is smaller than our BEG sample size, 
because the BEGs with each line S/N > 3 are 
plotted here. SF BEGs are not plotted in (b). 



2. OBSERVATION 
2.1. Target Selection 

We selected 59 non-passive (i.e. currently star- 
forming or ho sting an AGN) B EG targets based on 
the scheme of iLee et al.l ( 20081) . using the spectro- 
scopic sa mple of galaxies in the SPSS Data Release 
4 (DR4; lAdelman-McCarthv et ail 120061 ). We 
adopted the physical parameters of the galaxies 
from several value-added galaxy catalogs (VAGCs) 
drawn from the SDSS: photo metric parameter s 
from the SDSS pipeline (IStoughton et al.l [20 021. 



struct ural parameters e stima ted by IPark fc Choi 
(|2005l ) and IChoi et al.l (|2007l ). and spectroscopic 
parameters from Max-Planck-Institute for Astro- 
physics ( MPA) / Johns Hopkins University( JHU) 
VAGC (IKauffmann et all 120031 : iTremonti et~al 



2004HGallazzi et al.ll2006ll Firstly, we selected the 
early-type galaxies from the SDSS galaxies using 



their distribution in the color, color gradient and 
light- concentration parameter space (|Park fc Choi 



20051 ). Next, we divided the selected early- type 



ga laxies into REGs and BEGs, using the method 
of ILee et al.l (|2006l) based on the color distribu- 
tion of bright early-type galaxies as a function of 
redshift q The BEGs were classified using their 
flux ratio s between several spectral lines (e.g . BPT 



diagram; Baldwin et al.lll981t IKewlev et alj |2006) 
into passive, SF, Seyfert, and low ionization nu- 
clear emission region (LINER) galaxies, among 
which passive BEGs were excluded. Finally, we 
dropped BEGs with conspicuous sub-structures 
or disk components in visual checks. Among the 
BEGs selected in this scheme, we chose the target 



2 Here, BEGs are bluer than the Gaussian peak in the [g — r) 
color distribution of the bright ( 01 M r < —20) early-type 
galaxies, by more than 3ct, where a is the Gaussian spread. 
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Fig. 3. — The atlas images of the observed 36 
BEGs, retrieved from the SDSS. The bar in each 
image represents a scale of 2" or 5". 



objects as non-passive BEGs at 0.02 < z < 0.1 
with K s < 14.4 Q The magnitude limit was ap- 
plied to ensure the quality of the AKARI spec- 
troscopy. 

As a result, 59 BEGs were selected and pro- 
posed to be observed, but only 36 targets among 
them were successfully observed by the AKARI 
IS Fig. Q] shows the basic optical-NIR properties 
of the 36 BEGs in our sample: their magnitudes 
as a function of redshift, color - magnitude rela- 
tion and color - color relation, confirming that our 
sample BEGs are blue and relatively bright ob- 
jects. Fig. [5] displays the distribution of the BEGs 
in the line flux ratio diagrams, showing that the 
SF:Seyfert: LINER composition in our BEG sam- 
ple is balanced well (11:15:10). In Fig. El the atlas 
images of the BEGs retrieved from the SDSS are 
presented. The basic properties of the BEGs are 



3 The '20 mag arcsec 2 isophotal Ks fiducial ellipse aper- 
ture mag nitude' from the T wo Micron All Sky Survey 
(2MASS; iJarrett et ail l2000l) Extended Source Catalog 
(XSC). http://www.ipac.caltech.edu/2mass/ 

4 There is various unpredictability in the space telescope op- 
eration, which makes it difficult to carry out the observa- 
tions exactly as planned. Because of such unpredictability, 
the perfect success of proposed observations is not guaran- 
teed. 



summarized in Table [T] 

2.2. AKARI Spectroscopy 

The AKARI |£| is a Japanese infrared astronomy 
satellite, from the Institute of Space and Astro- 
nautical Science (ISAS) of the Japanese Aerospace 
Exploration Agency (JAXA). It was launched on 
21 February 2006 and ran out of its on-board 
coolant (Helium) supply on 26 August 2007 after 
its successful operation. The observation of our 
BEG targets was conducted during the AKARI 
open-time operation, the post-Helium phase of the 
AKARI observation, in which NIR (1.8-5.5 fim) 
imaging and spectroscopic observations are avail- 
able. 

The observation was conducted from Novem- 
ber 2008 to August 2009. The target BEGs were 
observed using the IRCZ4 Astronomical Observ- 
ing Template (AOT), which is the spectroscopic 
AOT of the IRC. Between the grism and prism, 
we selected the grism spectroscopy that has better 
spectroscopic resolution. The 1-er noise equivalent 
flux density of the IRCZ4 grism is about 0.2 mJy, 
and its spectroscopic resolution is R ~ 120. As 
the target positioning option, we used the point 
source aperture (Np), because the typical effec- 
tive radius of target BEGs is only several arc- 
seconds. Among the successfully-observed 36 tar- 
gets, 29 objects were observed twice, while seven 
targets only once. The net exposure time for a 
single pointing observation is about six minutes 
(|Onaka et aLll2007t limanishi et al.ll2008l ). 



3. DATA REDUCTION AND ANALY- 
SIS 



3.1. Basic Data Reduction 
Stacking 



and Median 



Unlike our expectation, two pointings per each 
object were not enough to secure sufficient signal- 
to-noise ratio (S/N) for most targets (a half of 
the spectra have 3.5 — 3.9/im continua with S/N 
< 6), although the spectra of a few objects are rel- 
atively good. Thus, we applied the median stack- 
ing analysis technique for high S/N. The steps of 
the stacking analysis are as follows: 



5 http: //www. ir.isas.jaxa.jp/ASTRO- 
F/Outreach/index_e.html 



Step.l We used the IRC Spectroscopy Toolkit for 
Phase 3 data Version 20090211 provided by 
the AKARI team, for the basic processes 
such as dark substraction, image flat field- 
ing, spectral flat fielding, image combining, 
source detection and extraction, sky sub- 
straction, wavelength calibration and flux 
calibration. 

Step. 2 Using the SDSS spectroscopic redshift, we 
converted the observer-frame wavelength 
into the rest-frame wavelength for each ob- 
ject. Since the rest-frame wavelength ranges 
are different between objects, we trimmed 
each spectrum with the rest-frame wave- 
length range of 2.5 - 4.5 /iin. 

Step. 3 All trimmed spectra were normalized for the 
median flux density at 3.5 < A < 3.9 /im to 
be identical. 

Step. 4 We stacked the normalized spectra by taking 
a median. We also produced several median- 
stacked spectra in bins of several properties 
such as optical spectral type, optical color, 
and [Oii]-based specific star formation rate. 

We selected the median stacking instead of 
mean stacking, because it removes unreal data 
points (hot or bad pixels) more effectively, whereas 
it improves data quality similarly to that the 
mean stacking does. Not only the total stacking 
(i.e. stacking all spectra into a single spectrum) 
but also subsample stacking was carried out to 
compare the NIR spectral features between sub- 
samples with different optical properties. We sim- 
ply estimated the uncertainty of the stacked spec- 
tra, by supposing that the spectral noise follows 
the Poisson error and that the S/N is accumulated 
during the stacking process. That is, 



(S/N) 8 = JJ2(S/N) 



(1) 



where (S/N) s is the S/N of the stacked spectrum 
and (S/N), is the S/N of the individual spectra. 



3.2. Spectral Fitting 

We used the power law continuum 



PAH emission (Gaussian) formula ([Risaliti et al. 



3.29 /mi 



2006) to fit the stacked spectra: 

(\-3.29fJ,m) 2 



i — i — i — I — i — i — i — r 

EW=17.8±0.6nm 
r=-2.56±0.01 




[fim] 



Fig. 4. — The stacked spectrum (noisy thick solid 
line) using all BEGs in our sample. The noisy thin 
solid lines show the S/N-accumulated noise and 
the smooth solid line is the continuum + PAH 
emission fit using Equation (j2]). The estimated 
PAH equivalent width (EW3.29) and the contin- 
uum slope (r) are denoted on the upper-right cor- 
ner. The dashed vertical line shows the PAH cen- 
tral wavelength (3.29 /^m), and the dotted vertical 
line shows the upper limit of the wavelength range 
for fitting (3.9 /im). 



where F\ is the flux density, A is the wavelength 
in unit of //m, T is the continuum slope, a is 
the Gaussian spread in unit of /xm, and A and 
B are amplitudes. The fitting was conducted in 
the wavelength range of 2.5 - 3.9 //m, because 
the Bra emission and CO absorption lines may 
affect the NIR spectra significantly at A > 4/im 



anect tne in in spectra significantly at A 
(jlmanishi et al.ll2008l: [Risaliti et al.ll2010h . 



F x = AX r + Be' 



(2) 



4. RESULTS 

Fig. @] displays the stacked spectrum of all 
BEGs in our sample. In this spectrum, a clear 
PAH emission line is found, the central wave- 
length of which is 3.29/xm. The continuum is 
overall fit well to the function in the fitting range 
(2.5 — 3.9/Ltm). The ice-cored dust absorption (A ~ 





Fig. 5. — The stacked spectra of BEGs in three 
sub-samples with different [Oii]-based specific star 
formation rates (SSFRs). (a) S/Nr OII ] < 3 (that is, 
SSFR[ OII ] is not available); and among the BEGs 
with S/N [on] > 3: (b) SSFR [on] < 6.3 and (c) 
SSFR [on] > 6.3, in unit of 1CT 12 yr" 1 . 



2.8 — 3.2^m) feature is marginally found, but the 
bare carbonaceous dust absorption (A ~ 3.4/xm) 
feature is hardly seen in Fig. [4] Around A ~ 
3.4/im, instead of absorption features, a marginal 
excess is seen, which we currently have no idea of. 
It is noted that the long- ward of 3.9 fim of the 
stacked spectrum shows values totally lower than 
the fit. The known line features at A > 4/zm is 
Bra emission (A = 4.05/im) and CO2 absorption 
(A = 4.26/im), but those features are not clearly 
identified in Fig. [4] The apparent 'dip' at A > 4/xm 
may be simply because the spectral fitting was car- 
ried out in the wavelength range of 2.5 — 3.9/iin. 

Figs. [5] - [7] show the stacked spectra of the 
BEGs in bins of several properties: [Oii]-based 
specific star formation rate (SSFR[ OII j), 01 (u — r) 
color, and optical spectral type. In the estima- 
tion of SSFR, the star formation rate (SFR) was 
estimated using the [On] emission line that is 
know n to be hardly contaminat ed by AGN emis- 
sion (JHol 120051 iKim et alj|2006| ). and the stellar 
mass was derived using the 2MASS K s magni- 



Fig. 6. — The stacked spectra of BEGs in 
three subsamples with different optical colors, (a) 
ai (tt-r) < 2.2, (b) 2.2 < 0A (u~r) < 2.5 and (c) 
0A (u-r) > 2.5. 



tude (see iLee et al.1 l2010al) . When considering 
the objects with sufficient [On] emission signal 
(S/N [oil] > 3; Fig. [5b and Fig.Efc), PAH emission 
is hardly found for low SSFRr OII ] BEGs, whereas 
the stacked spectum of the high SSFRr OII i BEGs 
shows a clear PAH emission feature. This confirms 
the fact that PAH emission reflects current SF. It 
is noted that 14 BEGs have low S/Nr OII i, which 
may be mainly due to the difficulty in measuring 
the [On](3727A) line using the SDSS spectroscopy 
(3800 - 9200A) at low redshift. 

Fig. [5] shows that the bluer BEGs have the 
stronger PAH emission feature, which indicates 
that the blue colors of BEGs reflect their SF activ- 
ity. This trend is also confirmed in Fig. [71 the PAH 
equivalent width (EW3.29) of SF BEGs (25.2 nm) 
is larger than those of Seyfert and LINER BEGs 
(17.2 nm and 13.4 nm, respectively). In Fig. [3 it is 
noted that the EW3.29 of Seyfert BEGs (17.2±0.8 
nm) is larger than that of LINER BEGs (13.4±1.4 
nm) by more than 2a. Seyferts s how stronger 
AGN activity than LINERs (e.g. iGroves et al 
120061 ) and the current concensus on the relation- 
ship between AGN activity and PAH/SF is that 
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Fig. 7. — The stacked spectra of BEGs in three 
sub-samples with different optical spectral types, 
(a) SF BEGs, (b) Seyfert BEGs and (c) LINER 
BEGs. 



AGN activi ty destroys PAH particles and sup- 
presses SF dVoidll992t lAntonuccio-Delogu fc Silk 



20081 iRaffertv et al.ll2008l) . Apparently, however 
our results seem to be contradictory to such pre- 
vious understanding (i.e. Seyfert BEGs are more 
active than LINER BEGs in both AGN and SF 
activity) . 

5. DISCUSSION 

5.1. The Identity of BEGs from the AKARI 

View 

BEGs arc known to have not only very young 
stars but also a considera bl e amount of ol d stars 
(IMenanteau et al.l l2001ai iFerreras et al.1 l2005t 



Lee et al.ll2010al) . In our results, however, it is 



not easy to check this previous knowledge inde- 
pendently, because the NIR continuum slope (r) 
depends on both stellar and dust contents (dust ei- 
ther in host galaxies or in AGN tori). Fig. [8] shows 
the variation of T as a function of age, metallicity, 
star formation history and dust attenuati on, using 
the models of iBruzual fc Chariot! (J2003I ): simple 
stellar population (SSP), exponentially-decreasing 
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Fig. 8. — Comparison of NIR continuum slope 
(r) between population synthesis models and the 
stacked spectra of BEGs: (a) the simple stel- 
lar population (SSP) models with solar metallic- 
ity (Z=0.02) and various dust attenuation values 
(solid, long-dashed, short-dashed and dotted lines 
for Ty =0,1,5 and 10, where Ty is the dust opti- 
cal depth in the V band) ; (b) the SSP models with 
Ty — and various metallicity (short-dashed, solid 
and long-dashed lines for Z=0.008, 0.02 and 0.05; 
and (c) the exponentially-decreasing star forma- 
tion rate model with the exponential time scale 
of 1 Gyr (EXP; long-dashed and dotted lines for 
Ty — and 10) and the model with constant 
star formation rate of 1 M /yr (CONST; solid 
and short-dashed lines for Ty — and 10), with 
Z=0.02. The shaded area shows the T range of the 
BEG stacked spectra and the horizontal solid line 
shows the T value of the all-BEG-stacked spec- 
trum. Dust att enuation was es t imate d using the 
simple model of ICharlot fc Falll (|2000l ). 



and constant star formation rate models. In Fig. [51 
it is found that T becomes less negative as star for- 
mation ended earlier, as age increases, as metal- 
licity increases and as dust increases. Due to 
this degeneracy, it is difficult to tell the domi- 
nant factor in determining T of BEGs. Even if we 
suppose that the NIR continuum of the BEGs is 
old-SSP-dominated from the previous knowledge 
(the age of most stars in BEGs is about 10 Gyr; 
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Fig. 9. — The internal extinction E(B—V)i versus 
the r and EW3.29 of the BEGs. Filled circles are 
the values from the subsample-stacked spectra and 
cross is the value from the all-BEGs-stacked spec- 
trum. The horizontal errorbars indicate the sam- 
ple inter-quartile range for the median E(B — V)j 
value in each subsample. 



e.g. iFerreraset all 1200a iLee et alll2010aD . the T 
of the BEGs can be reproduced using either the 
metal-rich SSP model without dust or the sig- 
nificantly dust-attenuated SSP model with solar 
metallicity. 

The effect of internal dust extinction is inde- 
pendently checked using the information extracted 
from the SDSS data; that is, usi ng the Balmer 
decrem ent. Based on the formula of lCalzetti et all 
( 19941 ). we derived the internal reddening E(B — 
V)i (median value for each subsample) and com- 
pared them with the T and EW3.29 of the stacked 
BEGs in Fig. M It is found that the NIR spec- 
tral parameters have correlations with E(B — V)%, 
in the sense that the T and EW3.29 increase with 
increasing E(B — V)i- The E(B — V)i values 
of the BEGs are not small, but seem to be not 
large enough to explain the T of the BEGs in 
Fig. [8] ( the t v corresponding to E(B - V) = 0.4 
is 1.14; |Pe Ruvter et alJlJOOa ). Thus, if we sup- 
pose moderate (not too strong) dust attenuation 
in the BEGs, Fig. [^indicates that the BEGs may 



Fig. 10. — Comparison between BEGs and other 
objects in the T - EW3.29 plane. The data for 
starburst ULIRGs (SB ULIRGs; open rectangles), 
AGN ULIRGs (filled rectangles), LINER ULIRGs 
(filled triangle s) and not-cla s sified ULIRGs (open 
stars) are from lRisaliti et all (|2010h . The filled cir- 
cles are BEGs stacked in bins of several properties 
and the cross shows the all-BEGs-stacked result. 



have stellar populations with slightly high metal- 
licity, although this estimation is strongly model- 
dependent. 

Meanwhile, BEGs are often thought to be the 
intermediate objects transforming from interact- 
ing galaxies or mergers to re d elliptical galaxies or 
bulges of late-type galax ies (|Lee et allbool |2008t 



Kannappan et all 120091 ) . Known to have simi- 



lar origins are Ultra-Luminous InfraRed Galax- 
ies (ULIRGs). ULIRGs are very luminous in 
the IR band due to vigorous starburst a nd /or 
AGN activity (jVeilleux et all Il995l Il999ah . and 
are thought to be formed by strong interaction o r 
merger of two disk galaxie s (| Sanders et al.lH988l) . 
Risaliti et all (|2006l 120101) estimated the T and 
EW3.29 of some ULIRGs, showing an interesting 
division between AGN-host ULIRGs and starburst 
ULIRGs. 

Fig. [10] shows the loci of our BEGs on the F 
- EW3.29 plane, compa r ed wi th the objects pre- 
sented in iRisaliti et all (|2010l ). In this figure, it 



is found that our BEGs have T and EW3.29 val- 
ues distinct from those of the ULIRGs, in the 
sense that the BEGs have largely negative T values 
(~ —2.5), while those of the ULIRGs are mostly 
larger than -2 (up to 6). The EW3.29 of the BEGs 
is smaller than that of the ULIRGs on average, 
showing that the SF in the BEGs is not as active 
as that in the ULIRGs. The F value is expected 
to be largely positive when the NIR continuum is 
dominated by dust emission. In our BEG sam- 
ple, however, even the stacked spectra using opti- 
cal AGN or optical SF BEGs has largely negative 
r. This indicates that the dust amount in those 
BEGs may not be as large as that in ULIRGs, 
as shown already using the internal dust redden- 
ing derived from Balmer decrement (the mean 
E(£ - V) value of ULIRGs is larger than 1.0; 
Veilleux et aHll999al lbh. 

In short, the Y and EW3.29 features of the 
BEGs are not easy to interpret, because they are 
affected by the combined effects of age, metal- 
licity, star formation history and dust attenua- 
tion. However, based on some previous knowledge, 
Balmer line information and a few assumptions, 
the BEGs are thought to be old-SSP-dominated 
metal-rich galaxies with moderate dust attenu- 
ation. The dust attenuation in the BEGs may 
originate from recent star formation or AGN 
activity. This interpretation is consistent with 
the previous understanding of stellar populations 
in BEGs (t hat is, mostly old stars + partiall y 
young stars; iFerreras et al.ll2005tlLee et al.l2010al) . 
There is a possibility that ULIRGs are the progen- 
i tors of BEGs in th e merging/interacting phase 
(jVeilleux et al.1 120021 ) , but any clear evidence of 
the close relationship between BEGs and ULIRGs 
is not yet found in our results based on the NIR 
spectroscopy. 

5.2. NIR Features and AGN Activity in 
BEGs 

As introduced in §1, NIR spectroscopy provides 
interesting clues to the internal evolutionary pro- 
cess of the BEGs, related to their SF or AGN 
activity. Those clues become more useful when 
combined with the optical properties of the BEGs. 
Fig. [TT] compares the stacked spectra in bins of 
several properties on the T - EW3.29 plane, giv- 
ing a summary of the key results in Fig. [4] — [7] 
The EW3.29 reflects the current SF of the BEGs 
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Fig. 11— The V versus EW3.29 of the BEGs 
stacked in bins of several properties: (a) all, (b) 
°- 1 (u — r) color, (c) SSFR[ OII ] and (d) optical spec- 
tral type. In every panel, the all-BEGs-stacked 
result is denoted as a cross. 



well (Fig. nib), but it is noted that the Seyfert 
BEGs have EW3.29 larger than that of the LINER 
BEGs. As mentioned at the end of SJH this result 
is not well explained simply by the current con- 
census that AGN activity destroys PAH particles. 

In addition to EW3.29, another key parame- 
ter is r. As shown in Fig. HI V is affected com- 
monly by age, metallicity and dust. However, the 
bluest BEGs with the largest (least negative) F 
in Fig. [TTb is not easily understood by the age or 
metallicity effects, because blue galaxies tend to 
be young or metal-poor, which will make Y more 
negative unlike the bluest BEGs. Thus, the Y dif- 
ference between the bluest BEGs and the other 
BEGs seems to be mainly due to the difference 
in their dust contents. It is also noted that the 
optical SF BEGs have small (largely negative) Y 
in Fig. Hid , which shows good consistency with 
the star formation history dependence of Y, in- 
dicating that those SF BEGs are not very dusty. 
That is, several AGNs with dusty tori in the bluest 
BEG subsample seem to mainly contribute to the 
largest (least negative) Y. 

In Fig. [TTb, the BEGs with large SSFR [on] 



have T intermediate between those of the bluest 
BEGs and optical SF BEGs. These are because 
the optical color, SSFR and optical spectral type 
of the BEGs do not tightly correlate: some non- 
SF BEGs have large SSFR [oil] or blue 01 (u - r) 
color. These results lead us again to the conclusion 
that the SF and the AGN activity in BEGs are 
not necessarily contradictory to each other. The 
coexistence of SF and AGN is not entirely new dis- 
covery (e.g. Tzanavaris fc Georgantopoulosll2007t 



Treister et al.1 120091 ) . but currently it is widely 
believed that AGN activity tends to destroy 



PAH particles a n d to t runcate SF ( Void 



Weinman n et al.l | 2006t iSchwanski et al 



1992; 



2007; 



Tortora et al.l 120091) . In our result, however 



Seyfert BEGs have larger EW3.29 than LINER 
BEGs, in spite of the fact that Seyferts show 
stronger AGN activity than LINERs. 

This resu lt may be exp lained by variou s ways . 



Imanishi et"aH |2008) 



Smith et al.l ([20071 ) and 
showed that PAHs can survive even if they are 
close to AGNs when the AGNs are not dust-free, 
which gives an answer to how the AGN activ- 
ity and PAH emission coexist. However, it docs 
not sufficiently explain why stronger AGNs have 
stronger PAH emission features. One possibility 
is that Seyfert BEGs and LINER BEGs may not 
be the separate branches in BEG evolution. Sat- 
isfying our results, one possibl e scenario of BEG 
evolution is as follows (see also lLee et al . 2010a): 

Ph.l SF BEG phase. Triggered by some mech- 
anisms (possib ly interaction or merger; 
Lee et al.ll2006l ). active (but not very dusty) 



SF occurs in the progenitor of a BEG, origi- 
nally consisting of old stars in the main. The 
NIR continuum is very similar to that of old 
passive galaxies, but the PAH emission is 
clearly detected. 

Ph. 2 Seyfert BEG phase. After the central black 
hole in a BEG increases its mass by gas ac- 
cretion, AGN activity starts. The SF in the 
BEG starts to be suppressed by the AGN 
feedback, but at this time, both AGN and 
SF activities coexist. The NIR continuum 
is contaminated by the AGN (r increases), 
and the PAH emission becomes weaker than 
that of SF BEGs. 

Ph.3 LINER BEG phase. The AGN feedback has 
truncated most SF by removing ambient gas 



in the BEG. Since this makes the gas ac- 
cretion into the central black hole stop also, 
AGN activity becomes weak as a natural re- 
sult. The NIR continuum is intermediate be- 
tween SF BEGs and Seyfert BEGs, and the 
PAH emission feature is weakest among the 
three phases. 

After these three phases, the BEGs may evolve 
into REGs or bulges of l ate-type galaxies via 
the phase of passive BE Gs (|Lee et al.ll20ol l2008t 
Kannappan et al.ll2009l ). The proposed scenario is 
consistent with the findings about the relationship 
b etween Seyferts and LINERs in several recent 
studies ( Schwansk i et al.ll2007L l2009t iHickox et al. 



20091 iLeeet alJl2010al ). 



6. SUMMARY 

We conducted an AKARI NIR spectroscopic 
survey of SDSS-selected BEGs. We secured the 
NIR spectra of 36 BEGs that are well balanced in 
their SF/Seyfert/LINER composition. We stacked 
the BEG spectra all and in bins of several proper- 
ties, such as color, SSFR and optical spectral type. 
This is the first presentation of the BEG NIR (2.5 
- 4.5 fim) spectra, from which we estimated the 
NIR continuum slope (r) and the equivalent width 
of 3.29 fim PAH emission (EW3.29). 

In the comparison between the estimated NIR 
spectral features of the BEGs and those of the 
model galaxies, the BEGs seem to be old-SSP- 
dominated metal-rich galaxies with moderate dust 
attenuation. The dust attenuation in the BEGs 
may originate from recent star formation or AGN 
activity and the BEGs have a clear feature of PAH 
emission, the evidence of current SF. BEGs show 
NIR features different from those of ULIRGs, from 
which we can not find any clear relationship be- 
tween BEGs and ULIRGs. 

We compared the stacked spectra of the BEGs 
in bins of several properties, confirming that the 
PAH emission is a good indicator of current SF. 
We found that Seyfert BEGs have more active SF 
than LINER BEGs, in spite of the fact that the 
AGN activity in Seyfert BEGs are stronger than 
that in LINER BEGs. One possible scenario sat- 
isfying both this result and the AGN feedback is 
that SF, Seyfert and LINER BEGs form an evo- 
lutionary sequence: SF — > Seyfert — > LINER. 
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Table 1 
Basic Properties of the BEGs 



SDSS ID 



Rcdshift 



°- 1 M r 

(1) 



K s 

(2) 



Optical 
Spectral Type 



l (u- 



SSFR[ OII ] 
(3) 



004803.10+160057.9 


0.050 


-21.681 


14.189 


LINER 


2.676 


0.000 


013244.86-102250.3 


0.032 


-20.876 


14.171 


Scyfcrt 


2.503 


— 


024847.62-000633.0 


0.025 


-20.392 


13.964 


SF 


2.287 


— 


031654.91-000231.1 


0.023 


-20.087 


14.293 


LINER 


2.532 


— 


080328.06+250609.5 


0.028 


-20.888 


13.886 


Seyfert 


1.563 


— 


082443.28+295923.5 


0.025 


-19.743 


13.553 


Seyfert 


1.998 


— 


083021.32+354902.4 


0.053 


-21.818 


14.129 


SF 


2.099 


11.989 


084437.81+325423.2 


0.032 


-20.808 


14.011 


Seyfert 


2.290 


25.838 


090455.34+335722.1 


0.044 


-21.566 


14.013 


LINER 


2.572 


0.306 


094814.71+574444.4 


0.030 


-20.382 


14.175 


Scyfcrt 


1.984 


— 


094835.61+390342.5 


0.071 


-22.515 


14.242 


SF 


2.526 


— 


100152.34+110419.7 


0.052 


-21.531 


14.291 


SF 


2.033 


14.902 


110559.02+585645.7 


0.048 


-21.997 


13.501 


Seyfert 


2.557 


1.804 


113003.57+571829.4 


0.036 


-20.708 


14.227 


Seyfert 


1.863 


37.973 


113629.10+542102.9 


0.063 


-21.939 


14.218 


Seyfert 


2.513 


9.908 


114721.61+522658.4 


0.049 


-21.858 


13.964 


Scyfcrt 


2.019 


55.042 


115110.60+451322.3 


0.036 


-21.113 


14.356 


LINER 


2.287 


2.687 


123250.91+503041.0 


0.063 


-22.204 


14.250 


LINER 


2.299 


6.698 


130238.74+452637.3 


0.024 


-19.878 


14.358 


LINER 


2.340 


— 


130457.90+534651.0 


0.029 


-20.859 


13.587 


LINER 


2.481 


0.014 


134456.31+024359.5 


0.077 


-22.287 


14.267 


SF 


2.299 


— 


141229.34+642236.4 


0.036 


-21.358 


13.765 


SF 


2.578 


— 


141239.32-020231.1 


0.074 


-22.175 


14.215 


Seyfert 


2.798 


— 


141718.34+464153.5 


0.038 


-20.839 


14.229 


Seyfert 


2.542 


1.952 


143222.69+565108.3 


0.043 


-21.713 


13.813 


SF 


1.925 


9.136 


144918.44+523312.2 


0.066 


-22.270 


14.113 


LINER 


2.457 


3.405 


145323.38+390413.6 


0.032 


-20.798 


14.086 


SF 


2.270 


— 


145455.44+453126.5 


0.037 


-20.924 


14.283 


Seyfert 


2.356 


0.000 


153344.45+391341.9 


0.040 


-21.176 


14.160 


SF 


2.463 


0.365 


153928.55+465859.7 


0.038 


-20.770 


14.220 


Seyfert 


2.509 


— 


154744.14+412408.2 


0.033 


-20.873 


13.954 


Seyfert 


2.299 


35.112 


163112.14+402912.1 


0.031 


-20.516 


14.279 


Seyfert 


2.476 


— 


164733.64+250707.8 


0.057 


-21.762 


14.328 


SF 


2.363 


6.926 


172611.22+591831.1 


0.027 


-20.241 


14.191 


SF 


2.128 


— 


210258.87+103300.6 


0.093 


-22.877 


14.289 


LINER 


2.583 


0.193 


225506.79+005839.9 


0.053 


-21.988 


13.905 


LINER 


2.359 


0.046 
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Note. — (1) The r-band absolute magnitude with K-correction (Petrosian magnitude), (2) 
the if s -band apparent magnitude (20 mag arcsec~ 2 isophotal K s fiducial ellipse aperture mag- 
nitude), and (3) the specific star formation rate in unit of 10~ 12 yr _1 (that is, the newly-forming 
stellar mass per unit stellar mass per year). For the objects with S/N[ OII ] < 3, their SSFRs were 
not estimated ( — ). 
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